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Abstract 


In this work, we present a new mathematical model for the spread of 
hepatitis C disease in two populations: human population and medical 
equipment population. Then, we apply the Jacobi wavelets method com- 
bined with the decoupling and quasi-linearization technique to solve this 
set of nonlinear differential equations for numerical simulation. 
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1 Introduction 


Viral hepatitis is a major health problem worldwide, comparable to that 
posed by other major communicable diseases, such as human immunode- 
ficiency virus (HIV), tuberculosis, malaria, or, more recently, coronavirus 
disease 2019 (COVID-19). In this work, we are interested in viral hepatitis 
C (HCV). 

Hepatitis C is an inflammation of the liver caused by the hepatitis C virus. 
The virus can cause both acute and chronic hepatitis. According to the fact 
sheet of the World Health Organization (WHO) updated on October 2017 
for hepatitis C, 71 million people have been estimated for chronic hepatitis C 
infection in the whole world, and approximately 399,000 people die each year 
from hepatitis C [28]. Until today, researchers could not develop a vaccine 
or effective treatment that heals hepatitis C at 100% [25]. 

The hepatitis C virus is transmitted by exposure to contaminated blood 
resulting from bringing the blood of an infected person into contact with 
that of a person likely to be contaminated directly (transfusion) or indirectly 
(equipment of contaminated injection for example). In 2016, WHO intro- 
duced global targets, for the care and management of HCV, a 90% reduction 
in new cases of chronic hepatitis C, a 65% reduction in hepatitis C deaths, 
and treatment of 80% of eligible people with chronic hepatitis C infections 
[30]. In Algeria, the president of the “SOS hepatitis association”, spoke in an 
interview about the need to draw up a national plan against viral hepatitis, 
which will aim to improve prevention, care, and the availability of drugs. 
He also mentioned that the prevention of viral hepatitis poses a problem 
in Algeria because there is no real prevention against these viral infections, 
especially at the dentist. It is obvious to know that the majority of con- 
taminations by these viruses are done during dental care. Therefore, raising 
awareness against viral hepatitis “B” and “C” is very important to detect 
these diseases, especially since they are silent. Indeed, better prevention re- 
quires better knowledge of the modes of transmission and the populations 
at risk in order to improve education and teach the appropriate protective 
measures. The last century has seen the emergence and rapid development 
of mathematical modeling, which plays an important role in assessing and 
anticipate the impact of Public Health programs. 

Over the last decade, a large number of mathematical models have been 
developed to simulate, analyze, and understand the dynamics of a population 
of hepatitis C. In a related research work, Martcheva and Castillo-Chavez 
[19] proposed a model to study the role of a chronic infectious stage on the 
dynamics of HCV over the long term. Incorporating the immune class in [10], 
in [32], the latency period was merged. In [4], the authors showed both the 
effect of processing and immigration. Another model describes the effect of 
isolating chronically infected people [15]. Several studies have been carried 
out in [11, 20, 23, 21, 22, 33] showing the impact of HCV treatment in drug 
users on the prevalence of the disease. The optimal control theory has been 
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used to understand the efforts made to prevent the spread of the disease by 
different measures and strategies [1, 31, 34]. 

Our aim in this article is to understand how hepatitis C disease can evolve, 
by highlighting the role of sterilization of infected material, modeled by or- 
dinary differential equation (ODE), unlike the model of Miller et al. [24], 
which targets the population of drug users. Therefore a single mode of con- 
tamination which plays the role of a vector of the disease and a single host. 
Our new model SIR-MI consists of taking into account other causes of con- 
tamination, such as dental equipment, toilet equipment, needles, tattooing, 
and piercing equipment in interaction with a mixed human population, and 
then we resolve this model. 

On the other hand, wavelet theory plays an important role in many areas of 
mathematics and applied sciences, for instance, signal analysis in medicine, 
image processing, signal processing, data compression, statistics, and numer- 
ical methods [7, 18]. In recent years, wavelets based on orthogonal polyno- 
mials have been used in many researches to solve different problems such as 
ODE, partial differential equations, fractional differential equations, optimal 
control, and variational calculus [2, 9, 8, 26, 27], and this is due to orthogo- 
nality property. We propose the Jacobi wavelets method with general indices 
(a, 8) in this work in order to obtain computational solutions. This method 
generalized other methods like Legendre wavelets and Chebyshev wavelets. 
The Jacobi wavelets method reduces an ODE to a system of algebraic equa- 
tions by using the operational matrix of the derivative of Jacobi wavelets. In 
our numerical simulations, we have found that using the operational matrix 
of derivative simplifies the implementation of the method compared to us- 
ing the operational matrix of integration [3]. Then, we apply the decoupling 
and quasi-linearization technique (DQLT) combined with the Jacobi wavelets 
method to solve the underlying problem. 

In this paper, we propose, in section 2, a mathematical model SIR-MI that 
describes the dynamics of a population of hepatitis C. Section 3 will concern 
the mathematical analysis of the proposed model. Section 4 is devoted to 
explaining the different steps that lead to the implementation of the Jacobi 
wavelet method combined with DQLT. In Section 5, we apply the proposed 
method to simulate the model SIR-MI. Finally, Section 6 presents our con- 
clusions. 


2 Model formulation 


In order to understand the effect of sterilization of the material on the trans- 
mission and dynamics of hepatitis C, we propose a mathematical model SIR- 
MI developed by Miller et al. [24] with five compartments. That is, let Ny 
be the total population of humans, which is subdivided into three subclasses: 
Sw (susceptible), I (infected), Ry (recovered). The total population of Nay 
material is divided into two subclasses: My (uninfected), My; (infected). 
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The graphical representation of the proposed model is shown in Figure 1. 
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Figure 1: The compartmental model diagram. 


The model SIR-MI is given by the following system of ODEs: 


ds BS (t)My(t) 
“EO Fa My eH SHO 
Un —_BSult) Milt) 


Fe) = Baty My ty 7 HO ~ (0 + oe) Tat), 


a ©) = \n(t) — pa Ru(t), 

dM aMy (t)I y(t) 

a = Am — Tul) + In) + 6Mr(t) — kMu(t), 
dM, aMyu (t) In (t) 


aE S )= Su) + In) 0M7(t) —kMr(t), 


Si7(0) = SH, 111 (0) = Ty, Rr (0) = Ry, Mu (0) = Msg,, Mr(0) = M,,, 
Sx, lH), Ru,,Mu,,M1, > 9, 


(1) 


with 
Am =kNu #0. 
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The parameters used in our model are defined in Table 1. 


Table 1: Definitions of parameters used in model. 


Parameters Description 


An birth rate of susceptible 
LH natural mortality rate of the human population 
B rate of interaction between susceptible humans and infected material 
T mortality rate due to the disease 
Xr disease cure rate 
AM birth rate of uninfected material 
Qa interaction rate between infected humans and uninfected material 
k rejection rate of infected or non-infected material 
0 sterilization rate of infected material 


3 Mathematical Analysis of the Model 


In this section and in the first moment, we will apply the theorem of Cauchy— 
Lipschitz to demonstrate the existence and the uniqueness of the solution of 
the system (1). Then we will study the behavior of this solution by going 
through the calculation of the points of equilibrium as well as the stability of 
these points. However, first, we note that the population of material Nj, (t) 
is constant. 


Indeed, 
Z 4 dNu _ dMu dM; 
Nw (t) = Mr(t) + Molt) = —*() = "+ - 
So 
dM dM 
Yt) + F(t) = Ay — k (M7 (t) + My (t)), 
dt dt 
dNyu 
t) = 0. 
a) =9 


We show that the human population is not constant. So we have 


Nu (t) = Su(t)+In(t)+ Rut) <= “ Bey = (t) 4 on (t) 4 on (t). 

Then 

TH (t) + SH (ty + SH (8) = Aw — war (Salt) + In (8) + Rur(t)) — rad), 
ant (t) = —TInlt) 
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3.1 Existence and uniqueness of a positive solution 


To study the existence and uniqueness of the solution of problem (1), we need 
to apply the Cauchy—Lipschitz theorem. 

Our model (1) is a system of nonlinear, autonomous first-order differential 
equations that can be written as the following Cauchy problem: 


X' (t) = F(X (t)),t € [0,7], 9 
‘x (0) = Xo, 2) 
with 
Su (t) A(X) 
Tx (t) fo (X (t)) 
x()=| Ra) | ana F(x) =| BX) |, 
Mo (t) fi (X (b) 
Mr (t) fs (X (t)) 
where Ceera? 
fu (0) = dn — HOMO Su), (3) 
fa(X (0) = HOMO _ sty) — (rt un) In, 4) 
fs (X (t)) = Ala (t) — we Ru (t), (5) 
AO) =u ae +M;(t) — k(t), (6) 
fs (X (0) = LMU OHO _ gn) — kar). (7) 


7 Sa(t) + In(t) 
We recall that the norm Norm(-) in the space of continuous functions 
from J to R® (denoted by C (I,R*)) is defined by 
Norm (F) = max||F ()ll2, 


with ||-||, is the usual Euclidean norm in R°. 
We are now able to state the following result. 


Theorem 1. The differential problem (1) admits a unique solution 
(Si (t), I(t), R(t), Mu(t), M;(t))” € R® for all t € [0,7]. 


Proof. To demonstrate that the Cauchy problem (1) admits a unique solu- 
tion, it suffices to show that the vector function F' of the equivalent problem 
(2), is Liptschizian. 

Let t € [(0,T], X1, X2 € R®. Then 
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(Xi (t)) — far (X2 (4))I, 
(X1 (t)) — fa (X2 (#))I, 
|F (X1(t)) — F (X2(¢))|| = max 4 | fs (41 (¢)) — fs (X2 (4), 
(Xi (t)) — fa (Xo (t))I, 
(X1 (t)) — fs (Xe (¢))| 


We assume that at any instant t € [0,7], the human population Ny(t) = 
Sa (t) + In(t) is between two real numbers strictly positive Nmin and NMnax- 

We will examine each of the components |f; (Xi (t)) — fi (X2 (t))|, ¢ = 
1,...,5. Therefore 


fa (X1 (¢)) — fr (Xe (£))| 


= [2S OM ayy (+ OSM 4 wa) 
5 |—Su, (t) Mr, (t) + Su, (t)Mr, (t)| + wa |-Sa, (@) + Sx, (tI. 
m(t) 


By adding and subtracting the term Sz, M7,, we have 


lf: (X1 (t)) — fr (Xe (4))| 
< aye (t) |-Mr, (t) + Mr,(t)| +8 


M1, (t) 
Nm(t) 


|—Si, (t) + Sx, (t)| 


M, 
Since Sz, < Nmax and = <1, then 
Nu 


Fi (Xi) = AXONS (Nina 8+ mir) I(t) — XalO. 


m(t) 


For (4) and following the same reasoning, we find 


fo (M1 (t)) — fe (X2 ()I 


< (yg Nna +B+A +r + mir) IXl6) ~ Xl 


The linearity of terms in (5) leads to 
| fs (X41 (t)) — fs (Xe (¢))| S (A+ wx) [X21 (4) — X2(4)I. 


From (6), we have 
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fs (X41 (t)) — fa XW) < nai |—Mu, (t)In, (t) + Mu, (t)In, (0)| 
HA 


+6|Mr, (t) — Mrz, (t)| + &|—Mu, (4) + Mu, (41. 


By adding and subtracting the term My, I7,, we have 


fa (X1 (t)) — fa (Xo (t))| S Nap Me (t) |-In, (1) + In, (6)| 
a 
Mate (t) |—Mu, (t) + Mu, (#)| 
+0 |.Mr, (t) — Mz, (¢)| + &|—Mou, (t) + Mv, ()|- 
Knowing My, < Nu, Na(t) > Nmin and x5 <1, we arrive at 


aNwu(t) 


min 


| fa (Xa (t)) — fa (X2 ()| < ( ta+0+k) || X(t) — X2(t)||. 


Finally, from (7) and following the previous steps, we have 


Ifs (X41 (O) — fs (X2(O)| S (stato 


min 


ba+o+k) |Xi(t) — X2(t)|| . 


Therefore, we have 


|F(X1 (0) — F(Xea(t))I| < C (|X) — Xe()II, 
with 


Nmax +B +A+T+ pH, A+ LH, 


C =max ( Nmax + 8+ HH, 


= —Z 
Nu (t) Nu(t) 
aN y(t) 


Nenin +a+o+k). 


3.2 Equilibrium points 


In this subsection, we will look for points of equilibrium Fy and EF, (Theorem 
2) and study their stabilities. We limit ourselves to the stability of the point 
Eo. The stability of the point EF, will be made numerically. 
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The basic reproduction rate Ro 


Understanding how an epidemic develops once it has appeared is crucial if 
we are to hope to control it. To do this, various models have been developed, 
which highlight the crucial role played by the Ro parameter, describing the 
average number of new infections due to a sick individual. As one can imag- 
ine, if this number is less than 1, then the epidemic will tend to die out, while 
it may persist or even spread to the entire population if Ro > 1 ([12]). 

We recall, for a given matrix A, that Sp(A) represents the spectrum of A 
and that the spectral radius of the matrix A, denoted p(A), is defined by 


p(A) = max {|A| ,A € Sp(A)}. 


The disease-free point is 


r 
(Sx =, In =0, Rx = 0, Mv = Nur, My = 0). 
HH 


We consider different infected populations of the model. That is, 


E(t) = oH _ Ny (t) — (7+ wa) Ind) 
and 
dM, Ls aMy (t)In(t) 
B= Soy tay 7 OMe) ~ bao). 


To be able to calculate Ro, we use two matrices F and V, where the 
matrix F represents the appearance of new infected; that is, what comes from 
other compartments and which enters the infected compartment following an 
infection, 


BSH 
F(Iy, Mr) = ee Nu 
Nu 


The matrix V represents all those who leave the compartments of the 
infected and those who come there for any other reason, 


V (In, Mr) = ae ag 


We have 
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A 
. N x +5 
_ at M HLH 
= aNubeH 0 
(A +r7+ Lun) AH 


The matrix —FV~! represents the next generation matrix. The basic 
reproduction rate is given by 


Ro = p(-FV~"). 


After calculating the eigenvalues of the matrix —FV~', we find 


= Ba = Ba 
At = 4/ Drrtunyerey ONd A2 = —4/ Ota RO): 


We then conclude 


_ Ba 
e= oa eT ®) 


The calculation of equilibrium points 


Theorem 2. The system (1) admits two equilibrium points Eo and EF}, for 
strictly positive parameters. They are given indeed as follows. 


1. If Ro < 1, then the point Eo exists and it is given by 


r 
Eo = (*4.0,0,.ar.0) 7 
LH 


Ro>1 and aB+apy > (k+8)(7 +A), 


then the endemic point EF exists and it is given by 


E, = ($4,104, Ry, Nu — M7, M7), 


with 
. An (T+ pH +A TZ 
St = ie 
Ii- Born — An (k4+0) (A474 “H) Nu 
7: aB (t+ pa +A)+ Ny (A+7 + pH) (ape — (Kk +0) (7 +A))’ 
Ba ag 
Mi; = Nu — M;, 
Mi = apuNuly 


((aun — (K+) (7 +A)) Ty + An (k +8) 
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Proof. The equilibriums of the system (1) are given by the solutions of the 
following system of algebraic equations: 


st M* 
pe eo 
St, MF 
pont AG; — (7 + pH) IT = 9, 
Nu 
Am — =U! + 6M¥ — kM¥ =0, 
e Si, + Ih 
Mil; k MF; =0. 


As the population of the material is constant Nyy = My + M7, then the 
system (9) is reduced to the following four equations: 


BSi,M} 


= pes 10 
AH Vu La Si = 0, (10) 
St,Mi 

Poet _ ue — (r+ um) Thy = 0, (11) 

M 
Ali, — HH Ry = 0, (12) 

7 Nu — Mt)It 

ou = Mir _ (4 9) Mt =0. (13) 


Si +15 
The sum of (10) and (11) gives 


a EE (14) 
LH 


From (12), it is clear that 
A 
Ri = — Ih. (15) 
HH 


To determine M7}, we replace (14) in (13) and obtain 


Mtg apaNmuly 
1 (Cop — (k +8) (7 +2)) Ti + An (k +8) 
Substituting (14) and (16) into (11), we find 


( Ba (An — (7 + HH + A) Ty) 
((apin —(k +0) (7 +A)) Ty + An (k + 9)) 


(16) 


(47+ a) Nar) Tj, = 0. 


(17) 
According to (17), we distinguish two cases: 
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Case I [7, = 0, we then find 


d 
Sy = —, Ry = Mp =0 et My = Ny. 
MH 


Hence the existence of the first equilibrium point Eo is as follows: 


mn 
Eo — (*4.0,0,.ar0) $ 
MH 


Case II 17, 4 0, we have 


Boa (Ag — (t+ wa +A) Ty) 
(Can — (k +6) (7 + A)) Te + An (k + 8)) 


(A+7+ un) Nm =0. 


We can easily write the last equation in the form 


BI? + Alt, =0 
with 


A=Ba(t+pat+A)+(A+7+ HH) Nu (apr — (k +8) (7 +4)), 
B=—-Ay (k+0)(A4+7+ pH) Nu + Badr. 


Hence, 


Born — An (Kk +0) (A +7+ Un) Nu 


Pa ce eae nana or 


Ti, = 


Then, the existence of the endemic point FE; is given by 


iy = (Si, Iz, Riz, Nu — Mz, M7), 


. . An—(T+HH +A) TR 
oe 


Fi Badu An (k +8) (A+7 + HH) Nu 
H aB(t +e +A)+(A4t7+ ue) Nu (ape — (k +0) (7 +A))’ 


Rig = —Tiy, 
LH 
Mi, = Nu — M;, 
Nu lt 
My = atHNutly 


(op — (h+ 8) (7 + A)) Tip + Aw (K+ 8)) 
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The positivity of the equilibrium points 


It is clear that the point Ep is positive (belong to the positive orthant) without 
condition. It then remains to show that the point EF is positive. This 
amounts to showing that S7,,[7,, Rj, Mj are positive. 


1. The positivity of S%, is according to (14). 


It is clear that the denominator of S7, is positive, so it suffices to study 
the positivity of the numerator 


AHQMEH + Aq (k +8) LH 
ab + apy —(k+0)(7 +A) 
a 


An —(T+HH +A) Ty = 


As the numerator S; is positive, then it remains to show that 5» is 
positive; that is, 


aB+apy > (k+0)(7T +A). 
2. The positivity of Ij, is according to (18). We let 
A 
Ty = =. 
2B 
We write A as a function of Ro defined in (8), 


R2-1 


Therefore A is positive if Ro > 1. 


We write B as a function of Ro, 


We note that B is positive if Rp > 1. 


3. The positivity of R}, is according to (15). We note that R}, is positive 
if Ij, is positive. 


4. The positivity of M7 is according to (16). We let 
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It is clear that M, is positive if Ij, is positive. Then, M7 is positive if 
Mp is positive. Indeed 


C 
Mz >0= — >0, 
C2 


with 
Cy = apy (A+7T+ bx) (a6 + apa — (k +6) (7+A)) 


+A (k+ 0) (abe + (7+ ba +A) aps) > 0, 
Cy = a8 (T+ he +A) 4+ (A4+74+ Ue) (Quen — (K+ 6) (7 +A)) > 0. 


We then deduce that Mg is positive if 


ab+apun > (k+0)(r+A), 


and therefore M7; is positive if Ij, is positive and if 


ab+apy > (k+0)(7T+A). 


3.3 Stability 


The stability of the equilibrium point [5] results from the stability of the 
Jacobian matrix of the system (10), (11), (12), (13) ( ie., its eigenvalues 
must be negative), which is given by 


J(SH, IH, Rx, Mz) 


Pr 0 0 _ oda 
“i . osm 
T 
_ Nu -Xv- (Tr + LH) 0 Nu 
0 » —LH 0 
—a(Nu = M,)In a(Nyu a M1)SH —aly (k i 0) 
(Su + In)” (Sux +I)" Sy +Iu : 


Theorem 3. It holds that Eo is locally asymptotically stable (the solutions 
must approach an equilibrium point under initial conditions close to the equi- 


librium point) if and only if 
Ro <i. 


Proof. The Jacobian matrix at point Ep is given by 
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BAn 
apie 0 _ 
ces 
O- Ate 0 = 
J (Eo) = ( ua) Nubw 
0 A HH 
N. 
( 2 wf =e 
AH 
-A 0 0 -B 
_{0-Co B 
7 0 D-—-A 0 
0 EF 0 -F 


We calculate the characteristic polynomial of J(Eo), 


det(J(Eo) — XI3) = —(X + ux)? (kA + kr + kun + OA 
+67 — 08 + On +(0+A4+7+ pH +k) X +X?) 
= —(X + px)” P(X). 


We have the first eigenvalues 


X, = X2=—pn <0, 


and 
P(X) =A+ BX +CX?, 


with 


A=kA+kr+kupn + 0\+ 67 + Onn — a8, 
B=(0+A4+7+ uA +k), 
C=1. 


Let us use Descartes’ rule [16] to show that the coefficients of the poly- 
nomial P do not change signs. 


It is clear that B and C are positive. It only remains to show that A is 


positive or equivalently 


1— Ro > 0. 


So, A is positive if Ro < 1. 


According to Descartes’ rule, the polynomial does not admit any positive 
root. Hence, the stability of the point Eo. 
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4 Method of resolution 


4.1 Jacobi wavelets 


The Jacobi polynomials Io") (a > —1, 6 > —1) are orthogonal polynomials 
on the interval [—1, 1] ([13, 29]) with the weight function 


w(a) = (L—a)*(1+2)", (19) 
where m is a positive integer, which represents the degree of the polyno- 


mial. These polynomials belong to the weight space L?, ({—1,1]). The Jacobi 
polynomials can be represented by the recursive formula given by 


Tre (a) 
_ (a+ 8+2m-—1) [a? — 6? 4 z(a+ B+ 2m)(a +B +2m—2)] 0,6) 
7 Im(a+B+2m—2)(a+B+m) m1 (2) 
(a+m—1)(@+m—-1)(a+84+2m) (4,8) 
m(a+B6+2m— 2)(a+8+m) Jin~2 (); (20) 
where 
TP) (x) =1, IO (x) = aaa ae ran e, (21) 


2 2 


As the Jacobi polynomials are orthogonal with respect to the weight func- 
tion w, then 


(He), HOO) =H bam, for alln,m EN, (22) 
cL? 


where 


2°+8+1P (a t+m+1)T(B+m+1) 
2n+1+atfB)mIT(a+8+m++1)’ 


noo fa = es 


on,m represents the Kronecker symbol, [ is the Euler gamma function, and 
(-,-)p2 denotes the inner product of L?,({[—1, 1). 


The Jacobi wavelets are defined by 


k+1 
252 7(4,8) (ok41 -1 
WP) (x) = Cr) Im”? (2**1¢—-In+1), BA<ae< sk (24) 
, ; otherwise, 
where k € N, n = 1,...,2* represents the number of decomposition levels, 


m = 0,1,...,M is the degree of the Jacobi polynomials (M € N*). The 
k+l 
coefficient Sue is for normality. 
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4.2 Decomposition in Jacobi wavelets basis 


Since the Jacobi wavelets family {uw ° a ) } | ok forms an orthonor- 


m>0 
mal basis in L?, ({0, 1]), we can express all functions f in L?, ((0, 1}) as a unique 
linear combination of elements of this basis: 


2" 0 
= Gat ©: (25) 
n=1m=0 
where Chm = (f,w (on8 ) . From the point of view of the numerical 
Li, ((0,1)) 


analysis, we take the truncated sum (its projection on finite space) 


> 3 Cn, Pe Sea (x). (26) 


n=1m=0 
Let 
T 
C= [c1,0; +++) €1,M,€2,0,-++,€2,M;-- +) C2k 05 dace , Cok | 5 
and let 
T 
a a,B a8 ap a,B a,B 
A °B) = Ee Deeg ne ’ so Vaal ee eee care mee | 
(27) 
We can find the following matrix notation: 
f(a) = CTU? (x). (28) 


In this case, the U9) are called the 2*(M +1) Jacobi wavelets vector and 
C is a 2*(M +1) vector. 


The operational matrix of derivative 


The derivative of the Jacobi wavelets vector U‘~:?) from (27) can be expressed 
by [17] 
dw (8) (x) 
dx > 
where D(“%) denotes the 2" (M +1) x 2" (M +1) operational matrix given 
by 


Di) ylo8) (x), 
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FOA 9 1... 0 
pera) O Pree 
: *s 0 
0 ... 0 F(@,8) 


FP) ig (M +1) x (M +1) matrix, where its (,7)th element is given by 


eae a Co (29) 


0, otherwise, 
in which pieeP ) and ek ) are defined from (23), and Os are given by 


ap) _ TG+8) (@G-YtatfS+)T(ath+)5) 
M1 OG +a + B) Pats) 


(30) 


3 yer (2(d+1)+a+/)T (a+d4+1) 


‘. T(6+d+2) 


d=j-1 


4.3 Description of the solution method 


In this subsection, we describe how to apply the Jacobi wavelets to solve 
ODEs. Then, we use the DQLT with Jacobi wavelets method to solve a set 
of nonlinear differential equations. In the end, we present the formula of 
errors calculation. 


Linear first order differential equation 


Consider the linear first order differential equation with initial condition 


f(a) + a(e)f(e) =9(2), 2 €]0,1), 
ae ou 


where fo is arbitrary constant. To solve the problem (31), we decompose 


f(a) in the Jacobi wavelets basis { Ges ” jr by estimating (28), 
1M 


n=1,..., 
m=0 


sey Ty 


f(x) =CTUO*) (2), (32) 
where C denotes the solution vector of the problem. Then, we have 


f'(z) = CT DOA GOP) (a). (33) 
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Now, by substituting (32)—(33) into problem (29), we get the following alge- 
braic system: 


CT DP) + a(w;)1)W (@;) = g(a), t=1,...,ne, (34) 


where J is the identity matrix and nc is the number of collocation points. 
We have to insert the initial condition 


fo = CTU*) (Q), (35) 


Equations (34) and (35) generate 2*(M +1) set of linear algebraic equations, 
which can easily be solved for the unknown C’ by using one of the method 
of resolution an algebraic system. Consequently, f(x) given in (32) will be 
easily calculated. 


Set of nonlinear differential equation 


To solve a set of nonlinear differential equations, we will use the DQLT to 
transform this problem by iterative steps into a set of decoupled and lin- 
earized differential equations, where each equation can be written as the 
problem (31). Then we use the Jacobi wavelets method described in the pre- 
vious subsection. Let us consider a set of p nonlinear differential equations. 
This iterative technique can be defined by 


Given initial profile f(° eae arere! f°, 
(FC) + a (a) FO = a (2, £0,199... A), 
oe I+ l U l 
(a) + aaa) A = 90 (2, HPP, FP Be), (36) 


141 
(ey faye Slade seein) 


where te and Fe are the approximations of the solution f; at the current 
and the precedent iteration, respectively. At each iteration, we apply the 
Jacobi wavelets method to solve p linear differential equation. Then, for 
(1 + 1)th iteration, we can calculate the decoupling error using the following 
formula: 


Epaur = max (||fi — fr lla, Ifa — fa" lass Ife — fot lla). (37) 


The procedure is terminated when the error of decoupling is sufficiently small. 
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Error estimation 


Since the ODEs solutions are only known at collocation points, the most 
appropriate norm is the euclidean norm if the exact solution is given. The 
accuracy of the proposed method is estimated by 


error = ||f(@) — fea(2)llp = » If (a4) — fea(as))’, (38) 


where fe, is the analytic solution, f is the approximate solution, and nc the 
number of collocation points. 


5 Numerical simulations of model SIR-MI 


In this section, we will study the stability of the point Ey, numerically. Then, 
we simulate our model to see the importance of studying the effect of the 
sterilization parameter infected material and management on the evolution 
of the human population. We apply the Jacobi wavelets with DQLT, which 
makes it possible to numerically evaluate the solutions of the ODEs and to 
build their graphs. We conclude our section with a discussion of the results 
obtained. 


5.1 The study of the stability of the second equilibrium 
point Fy 


The following table gives us the biological parameters that verify the condi- 
tions of existence and stability of the second point of equilibrium F): 


Table 2: The parameters verifying the stability of F1. 


Equilibrium point By, = (1622, 9920, 1668, 18172, 11828) 
Parameter AH B 4 Xr AM T Qa k 6 
Value 230 0.23 0.05116 0.086 1500 0.011 06 0.05 0.3 


1. For the conditions of existence, we have 
(Ro = 1.6313) > 1, (39) 


and 
(a8 + ay = 0.1687) > ((k + 8)(7 + A) = 0.0339). (40) 
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Hence, we have the existence of the point F}. 


2. For stability, the Jacobian matrix of system (34) at point Fy, after 
substitution of the parameters given in Table 2 is given by 


—0.1418 0 0 —0.0124 
0.0907 —0.1482 0 0.0124 
0 0.0860 —0.0512 0 
—1.5836 2.5875 0 —0.5778 


J (Ey) = 


The eigenvalues of J(£,) are 


vp, = —0.0512, 
Up2 = —0.6842, 
vp3 = —0.0918 + 0.05152, 
up4 = —0.0918 — 0.0515i. 


The eigenvalues of J(£,) have a negative real part, hence, the asymp- 
totic stability of the second equilibrium point F. 


Figure 2: The convergence of the system toward EF). 


Figure 2: The convergence of the system toward EF}. 


We note that the solutions obtained in Figure 2. All converge towards 
the equilibrium point E, when t > +00. 

Figure 3 shows that the five subpopulations converge after a fairly large 
time to the second equilibrium point F}. 


In what follows, we carried out simulation experiments with the parame- 
ters illustrated by Table 3. 
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Human subpopulations 


0 2 40 60 80 100 120 140 0 20 40 60 80 100 120 140 
Time by years Time by years 


(A) (B) 


Figure 3: Evolution of human and material subpopulations: (A) represents the human 
subpopulations, (B) represents the material subpopulations. 


Table 3: Variations in estimated values of biological data. 


The time 40 years 
L.C Sy In Ru My M, 
Value 60000 30000 20000 10000 20000 
Parameter An B Ly Xr Am T Qa ia 7) 
Value 230 = -0.072-0.6 0.01-0.05116 0.006-0.235 1500 0.011 0.1-0.6 0.05 0.17-0.4 


5.2 The impact of equipment sterilization on disease 
progression 


Infected humans 


2 30 36 40 


i 5 10 6 mo 8 0 6 40 a 6 10 15 20 
Time by years Time by years 


Figure 4: Evolution of human subpopulations for different values of 6. 


For different values of 6 = 0.2,0.3,0.4, we see, in Figure 4, the positive 
effect played by the sterilization parameter to reduce the number of infec- 
tions. This shows that better compliance with universal hygiene rules and 
recommendations for disinfection of nondisposable medical equipment and 
the development of equipment for use single should allow in the long term a 
quasi-disappearance of infections. 
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5.3 The impact of the transition rate from Iy, to Ry 


Infected humans 
Resistant humans, 


i i === 2-001 


Peni 9=0.015 
=r=9z0.02 
0 5 10 16 20 25 30 Ea 40 


Figure 5: Evolution of human subpopulations for different values of X. 


For different values of A = 0.01, 0.015, 0.02, the curves obtained in Figure 
5 have made it possible to understand the important role of good care for 
infected people. Being infected with HCV does not protect against the risk of 
a new infection, which could worsen the medical situation. The development 
of a better therapeutic strategy can significantly improve the quality of life 
of people infected with hepatitis C. 


6 Discussions and Conclusion 


World Health Organization recommends that countries develop national 
strategies to reduce the burden of disease associated with hepatitis C ham- 
pered by weak or lacking national surveillance systems and unreliable esti- 
mates of the burden of hepatitis C morbidity. 


In this work, we described and mathematically analyzed the dynamics of 
hepatitis C. The different numerical simulations were presented to see the 
behavior of the model at infinity, and the results obtained showed that the 
trends related to the prevention and management of infection considerably 
influence the subpopulations. We also applied the Jacobi wavelets method 
associated with the DQLT to obtain a numerical solution, which gave a very 
satisfactory results. 


Due to the lack of data, our model has not been validated for the case of 
Algeria. Nevertheless the results of this modest work constitute the bases of 
work to be continued and improved for a much more in-depth study. 
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